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Abstract The classical 1/Va and PC methods of con- 
structing binned luminosity functions (LFs) are revis- 
ited and compared by graphical analysis. Using both 
theoretical analysis and illustration with an example, 
we show why the two methods give different results 
for the bins which are crossed by the flux limit curves 
L = Liim{z). Based on a combined sample simulated 
by a Monte Carlo method, the estimate of two meth- 
ods are compared with the input model LFs. The two 
methods give identical and ideal estimate for the high 
luminosity points of each redshift interval. However, 
for the low luminosity bins of all the redshift inter- 
vals both methods give smaller estimate than the in- 
put model. We conclude that once the LP is evolving 
with redshift, the classical binned methods will unlikely 
give an ideal est i mate over the total luminosity range. 
Page fc Carreral (l200nh noticed that for objects close to 
the flux limit ^i/y^ nearly always to be too small. We 
believe this is due to the arbitrary choosing of redshift 
and luminosity intervals. Because 4'i/Va is more sen- 
sitive to how the binning are chosen than (f>pc- We 
suggest a new binning method, which can improve the 
LPs produced by the 1/14 method significantly, and 
also improve the LPs produced by the PC methods. 
Our simulations show that after adopting this new bin- 
ning, both the 1/Va and PC methods have comparable 
results. 
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1 Introduction 

Prom shortly after the first quasar was found until the 
present, considerable effort has been spent in obtaining 
samples to investigate their luminosity distribution as 
a function of redshift, known as the luminosity function 
(LP). The LP is very important because its shape and 
evolution provide constraints on the nature of activity 
and the cosmic evolution of quasars/active galactic nu- 
clei (AGNs). Up to now many statistical approaches 
have been proposed to investigate the LPs. These in- 
clude parametric techniques which assume analytical 
form for the LFs, and non-parame tric methods w hich 
usually need binning the data (see |Johnstonll201ll for 
an overall review) . 

Among the non-parametric methods, the esti- 
Avni fc Bahcailll980l:lEaleslll993tlEllis et al 



mator (see 



1996h is the most classical binned method and is par 
ticularly prevalent for its simplicity. Although more 
than four decades have passed since i ts original ver - 
sion (i.e., the famous 1/Vmax estimator. ISchmidtill968( ) 



was presented, the 1/Va method is not outdated and 
continues to be widely used in the literature (see 



Civano et al. 2011; Mao et al. 


2012; Marchesini et al. 


2OI2I: McAbine & Jarvis 2011 


; Padovani et al. 


2011; 


Patel et al.ll2012: Hiroi et al.ll2012: Yuan & Wane 


2012; 


Marcha & Caccianiea 20131 for latest use). On the 



other hand, authors have pointed out that the 1/Va 
method introduces a significant error for objects close to 
the fl u x limit (e.g.. Page fc Carreral2000t Cara fc Lister! 



2008[ ). IPage fc Carrera ( 2000f ) presented an improved 



method (hereafter the PC method) of constructing the 
binned LP. According to the result based on a Monte 
Carlo simulation, the authors believed their method is 
superior to the 1/Va method in many aspects. How- 
ever, their simulation was performed using a luminosity 
function which is unchanging with redshift (i.e. no evo- 
lution) and a single flux limit. The assumption of no 
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evolution is too particular to be able to represent more 
general situations. Furthermore, they focused on the 
situation of single sample (a single flux limit), and the 
discussion on multiple samples was not sufficient. In 
this paper, we revisit the l/Va and PC methods to find 
the reason of systematic error using graphical analysis. 
The situation when multiple samples are combined to 
obtain a LF is particularly discussed. By the way, in re- 
cent years some more s ophisticated and rigorous uieth- 



Survey B 



Survey D 



ods have emerged (e.g . . lSchafeij|2007l:lKellv et al.ll2008 



Christlein et al.l l2009t iTakeuchil 120101 : iJohnston et al, 
20121) . Nevertheless, it needs time for the new methods 
to be recognized and in widespread use. During this 
time, specifying deficiencies of the old methods could 
be helpful. 

Throughout the paper, we adopt a Lambda Cold 
Dark Matter cosmology with the parameters flm = 
0.27, nA = 0.73, and Hq = 71 km s^^ Mpc-\ 



2 Methods 
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Fig. 1 Sample B and D. All data points are simulated. 
Sample D is assumed to be deeper in all frequencies than 
Sample B. Let 5*;^^ and Sj^,-^ denote their flux limits re- 
spectively, then Si^rn < S^rn ■ The red solid circles represent 
objects with fluxes between and Sj^^, while the black 
solid squares represent those with fluxes above S^^. 



The differential LF (j){L,z) is defined as the number 
density of target sources per unit comoving volume per 
unit luminosity interval, i.e. 



dVdL 



iL,z). 



(1) 



More often it is defined in terms of log L. In this paper 
wc do not differentiate L and log L strictly. 

2.1 The l/Va method 

Th e l/Vg method o r iginat es from the celebrated paper 
by Avni fc Bahcall ( 1980l) . which generalized V/Vmax 
( Schmidt 19681 ) for multiple samples. Here we consider 
two flux-limited samples observed at the same frequen- 
cies. For simplicity, we assume the two samples to be 
not overlapping in survey regions, as shown in Fig. 1. 
Sample D (D denotes deep) is assumed to be deeper in 
all the sample frequencies than Sample B (B denotes 
bright). 

spectively, then 5;^^ 

If N objects appear in the interval ALAz {Li < L < 
L2, zi < z < Z2) around the bin center (L, z), the LF is 
estimated as 



Let Sl^„^ and 5';^,„ denote their flux limits re- 

lim liin' 



1 ^ 1 
— V— . 

1=1 ° 



(2) 



According to the scenario of lAvni fc Bahcalll |l980l ). 
the sample B and D can be combined into a "coherent 
sample" . It can be regarded as a single sample, in which 



each object is allowed to be distributed anywhere within 
the total volume. Thus the available volume Va for any 
object i in sample B and D can be calculated as 



S=B,D 



■top fjy 



(3) 



where Q.s{zi, z, 

dians of the S th survey. We have ^s{zi < z, 



iS 
max 



) is the effective survey area in stera- 

) = Oc 

max/ ' 



and fls{zi > zl^^^) = 0. 



) = "s 
ri_B and flo are the solid 



angles subtended by B and D samples on the sky re- 
spectively, zlgp is deflned as 



4op = min[22,2'L 



mm[z2,z{L^,Sii„,)], 



(4) 



However, one must keep in mind that the sources 
with fluxes between S^^^ and Sj^^ (represented by red 
solid circles in Fig. 1, hereafter red sources) can never 
appear in sample B. In the L — z plane (shown in 
Fig. 2), these sources are located between the curves 
L = Lf^^{z) and L = Lj^mi^)^ which represent the flux 
limits of survey B and D respectively. On the other 
hand the sources with fluxes above (represented by 
black solid squares in Fig. 1, hereafter black sources) 
may appear both in sample B and D. The above discus- 
sion helps us to distinguish the actual surveyed regions 
corresponding to red and black sources. 

2.2 The PC method 

The key point of the l/Va method is that it takes into 
account the contribution of object i to the number den- 
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sity o f the bin ALAz as 1/{ALV^). iPage &: Carrera 
jresented an i mproved method of constructing 
the binne d LF (also see [ Tzanavaris fc Georgantopoulos 



[2OO8; Yencho et alll2009t) . That is. the LF at the cen- 



ter of a bin with a luminosity interval Li and L2 and a 
redshift interval zi and Z2 can be estimated as 



N 



J Li <J z\ az 



(5) 



where N is the number of sources detected within the 
bin, and the double integral corresponds to the shaded 
area of bin 1 in Fig. 2. The key point of the PC method 
is to consider the actual surveyed region of a bin as a 
four-dimensional polyhedron in the volume-luminosity 
space, and the 4>pc for this bin is the four- number- 
density. 



> 




Fig. 2 L-z plane for the simulated sources from sample B 
and D. Only a part of the sources are plotted. The black 
solid line L = Lfi^{z) is the flux limit curve of sample B, 
and similarly for L = Lf^^{z). The shaded regions marked 
by 1, 2 , 2 and 3 represent the surveyed regions of bin 
1, bin 2 and bin 3. A few example sources with difi'erent 
locations and status, i.e. i, j, k, m, and n, are shown in red 
and black dots. The red dotted lines are auxiliary lines to 
illustrate the values of ziR^ and ziP„„. 



2.3 Comparison of the two methods 

Here we take into account three situations represented 
by bin 1, bin 2, bin 3 in Fig. 2. For them, the shaded 
regions are the regions of the volume-luminosity plane 
in the interval ALAz that has been surveyed (also see 
Page fc CarreraluOOOl ). and are marked by 1, 2 ,2 and 
3. For bin 2, because either red or black sources appear 



k m n ^1 



Fig. 3 Illustration of the surveyed regions and the area 
equal to [ALVa]- (a) illustrates the surveyed regions of 
bin 1, bin 2 and bin 3. (b) illustrates the area equal to 
[ALFi], [ALV;*=1, [ALK™] and [ALKT] of the ex- 

ample sources i,j,k,m and n respectively. 

in it, 2 and 2 are used to represent different surveyed 
regions. In the volume-luminosity space, these shaded 
regions are four-dimensional polyhedrons illustrated in 
Fig. 3(a). Their four- volumes can be calculated as 



Ai = Qi 



A,, 



D 



L2+AL 



^2" = i^D + f^s) 



A3 = {nD + ^B) 



L2 + AL 



L2+3AL 



L2+2AL 



dzdL. 



——dzdL. 

dz 
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-^dzdL. 

dz 



——dzdL. 

dz 



(6) 



(7) 



(8) 



(9) 



Therefore, according to section 2.2, the estimated LF 
by the PC method (j)pc for bin 1, bin 2 and bin 3 are 
calculated as Ni/Ai, iV2/(A2' +A2" ) and ^^3/^3, where 
A^i, N2 and A''3 arc the numbers of sources detected 
within bin 1, bin 2 and bin 3, and Ai, A2' , A2" and A3 
are defined by Eq.(5)-(8). It is emphasized that N2 is 
the sum of red and black sources detected within bin 2. 

In the use of the 1/Va method, the available volume 
Va of every source should be calculated. The calculation 
of Va for a source depends on its location in a bin. In 
Fig. 2, a few example sources with different locations 
and status are labeled as i, j, fc, m, and n respectively. 

1. bin 1: Only red sources can appear in the bin. For a 
source i (represented by the red solid circle labeled 
'i' in Fig. 2) in bin 1, as z^^^ < zi, we have 



D 



- dV 
dz 



dz. 



(10) 
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similar to source k with 
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(b) 
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Fig. 4 Available volume Va for an object i within bin 1. 
The area equal to [ALVJ] is indicated by the shaded region 
in (a). It is noticed that this area depend on the position of 
object i within bin 1. (b) and (c) show how this area changes 
dramatically in two extreme situations when an object i is 
located at the left and right margin respectively within the 
bin. 

The area equal to [ALV^*] for an object i (repre- 
sented by a red spot) is shown in Fig. 4(a). This 
area is clearly not the same as that of the surveyed 
region (the shaded region of bin 1, marked as 1 in 
Fig. 2). The deviation depends on the position of 
an object i in bin 1. If the object locates at the 
left margin of bin 1 (sec Fig. 4(b)), the area equal 
to [ALI/j'] (the shaded region in Fig. 4(b)) will be 
much smaller than that of the surveyed region (the 
shaded region of bin 1, in Fig. 2), and vice versa (see 
Fig. 4(c)). Hence 4>i/Va ^"^^ 4'pc will give different 
estimates of 4> for this bin. 
2. bin 2: Both red and black sources appear in the bin. 
However they belong to different surveyed regions 
(represented by 2 and 2 respectively) . For a source 
j (represented by the red solid circle labeled 'j' in 



Fig. 2) in bin 2, when 
have 



< zi and Z2 < z: 



dV 
dz 



dz. 



(11) 



For a source k (represented by the red solid circle 
labeled 'k' in Fig. 2) in bin 2, when zi < z^^ < z^ 
and Z2 < z^dxi '^^ have 



kB 



yK 



——dz - 
dz 



(12) 



Although source m (represented by the black square 
labeled 'm' in Fig. 2) is a black one, its situation is 



fit 



—dz 
dz 



n 



D 
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dV , 
——dz. 

dz 



(13) 



The total surveyed region of bin 2 is the sum of 
2 and 2". The area equal to [ALV^] (see Fig. 3(b), 
j specifically) is clearly smaller than that of the sur- 
veyed region (Fig. 3(a), 2 -1-2 ). Consequently, the 
density contribution of a source like j to bin 2 is pos- 
itive, and possibly leads jy^ to give a exaggerated 
estimate for bin 2. The area equal to [ALV^*^] as well 
as [ALVJ"] arc approximations of the surveyed re- 
gion of bin2. Overall, (/'i/v„ possibly gives a higher 
estimate of than 4>pc for this bin. 
bin 3: Only black sources appear in the bin. For a 
source n (represented by the black square labeled 'n' 
in Fig. 2) in bin 3, as Z2 < z'^ax < ^"^axi we have 



dV 

dz 



dz + fli 



dV_ 

dz 



dz. 



(14) 



The area equal to [ALVJ'] for an object n is the same 
as that of the surveyed region and hence 4>i/Va ^-^id 
(j)pc will give the same estimate of (f> for this bin. 

To sum up, the two methods give different results 
for the bins which are crossed by the flux limit curves 
L = Lum{z). 
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Fig. 5 L — zph 

established bv lYuan fc Wand (|2012l '). The sample consists of 
four surveys, each represented by different symbols. When 
these data points are binned to use the 1/Va or PC methods, 
much more bins are crossed by the flux limit curves L = 
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Fig. 6 This is an example bin of combined sample from 
more than two samples. It is divided into multi-regions by 
the flux limit curves, (a) and (b) illustrate the surveyed 
regions and the area equal to respectively. 



2.4 Multiple Samples 

Practically many samples with different depth, extent 
and degrees of completeness are combined to obtain a 
LF that spans a wide range of luminosities and red- 
shifts. When the sources from different samples are 
plotted on the redshift-luminosity plane and binned 
to use the 1/Va or PC methods, much more bins are 
crossed by one or more flux limit curves. Some bins are 
divided into multiple regions by the flux limit curves 
L = Lum{z) (see Fig. 5 for a example). In this case, 
the area equal to [ALV^*] for an object i is a gross ap- 
proximation to that of the surveyed region (see Fig. 6). 
4>pc and 01/ Va '^iU give different estimate of (j) for these 
bins. 

2.5 Applying the methods to real data 

In this section we apply the 1/14 an d PC methods 
to the combined sample established by lYuan fc Wang 
(j2012l ). Fig. 5 shows the L — z plane of the 
combined sampl e. The four sub-samples, MRCl 



( McCarthy et al.l Il996l) MS 4 jBurgess fc Hunstead 



20061) ■ BRL |Best et al.ifl999t) and 3CR R jLaing et al 



19831) are described in the paper of lYuan fc Wane 
20121 ). For comparison, the radio luminosity function 
(RLF) at 408 MHz estimated by the two methods are 
plotted together for all the redshift bins in Fig. 7. It 
is not surprising that the two methods give the same 
results at the bright end of the RLFs, corresponding 
to the situation of bin 3 discussed in section 2.3. It 
is clear that 4>i/Va a-^d (ppc give different estimates at 
the faint end and middle of the RLFs. trends 
to give a smaller estimate than 0pc at the faint end, 
while it gives a larger estimate in the middle of the 
RLFs. Especially for the high redshift bins, 4>i/Va gives 
a significantly larger estimate than (j^pc in the middle 
of the RLFs. Because the situation like that for source 
j in bin 2 (discussed in section 2.3) is more prevalent 
for high redshift bins. 



2.6 Monte Carlo simulation 

In this section we use a combined sample simulated by 
a Monte Carlo method to further compare the 1/14 
and PC methods. The simulation is performed us- 
ing a double-power-law model radio luminosity func- 
tion which is changing with redshift in the form of 
(1 -I- z)'^ . The Monte Carlo simulation produces four 
flux-limited samples containing more than 1 000 000 
sources. The flux limits of the four simulated samples 
are 0.6 Jy, 1.2 Jy, 2.2 Jy, and 3.0 Jy respectively, and 
the solid angles subtended by them are 0.6 Sr, 1.2 Sr, 
2.2 Sr and 3.0 Sr respectively. Without loss reality, a 
random spectral index following a normal distribution 
with an average of 0.7 is arranged for all the simulated 
sources. The four simulated samples are then combined 
into a "coherent sample" and binned LFs are produced 
for it in a range of redshift intervals using both meth- 
ods. These are shown in Fig. 8, on the left and 
4>pc on the right. The input model LFs take values 
of z = 0.12, 0.35, 0.75, 1.25, 1.75, 2.25 and are shown as 
dashed lines. 

As expected from section 2.3, the two methods give 
identical and ideal estimate for the high luminosity 
points of each redshift interval. However, for the low 
luminosity bins of all the redshift intervals both meth- 
ods give smaller estimate than the input model. This 
situation is especially serious for the 1/14 method. For 
the lowest luminosity bins of the 1.0 < z < 1.5, 1.5 < 
z < 2.0, 2.0 < z < 2.5 redshift intervals 4)i/Vc is signifi- 
cantly smaller than the input model. 

T he above results are no t in completely agreement 
with IPage fc Carreral |2000l ). Their simulation showed 
that (ppc is always a good representation of the input 
model over the total luminosity range. This is because 
their simulation was performed using a LF which is un- 
changing with redshift (i.e. no evolution). Once the 
LF is evolving with redshift, the classical binned meth- 
ods will unlikely give an ideal estimate over the total 
luminosity range. In Fig. 9 we show how this happens. 

In Fig. 9 (a), binl,bin2,bin3 are three bins located 
at different luminosity. Each bin is further divided 
into nine sub-bins. We assume each sub-bin is small 
enough and the four-volume involved is unit. The num- 
ber noted in each sub-bin represents the source number. 
Thus the variation of number across sub-bins indicates 
how the density changes along luminosity and redshift. 
According to the PC method, the density (i.e., LF) at 
the center of bin2 and bin3 is estimated as 10 and 4 re- 
spectively. This is in agreement with the input model 
density at Zmid (see Fig. 9 (b)). Binl is divided into two 
parts by the flux limit curve L = Lumiz). Owing to the 
positive evolution of density along redshift and decline 
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27.5 28.0 28.5 29.0 29.5 28.0 28.5 29.0 29.5 28.5 29.0 



log.„{L.n./WHz') 

^10 ^ 408 f 

Fig. 7 Comparison of the RLFs at 408 MHz estimated by the binned 1/Va and PC methods are plotted together for aU 
redshift bins. The two methods give same results at the bright end of the LFs, but they give different estimates at the faint 
end and middle of the RLFs. 
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Fig. 8 Binned luminosity functions of simulated samples of objects using (left) 4>i/Va ™d (right) i^pc- The input model 
LFs take values of z = 0.12, 0.35, 0.75, 1.25, 1.75, 2.25 and are shown as dashed lines. 
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Fig. 9 This figure show that as long as the LF is evolv- 
ing with redshift, the classical binned methods will unlikely 
give an ideal estimate over the total luminosity range, (a): 
binl,bin2,bin3 are three bins located at different luminosity. 
Each bin is further divided into nine sub-bins. Each sub-bin 
is assumed to be small enough and the four- volume involved 
is unit. The number noted in each sub- bin represents the 
source number. Thus the variation of number across sub- 
bins indicates how the density changes along luminosity and 
redshift. (b): The model density (LF) at Zmid is shown as 
black dashed lines. The estimated density of the three bins 
by the PC method are shown as black dots. 



of density along luminosity, the left part of binl proba- 
bly contains more sourees than the right part. However, 
the sources in left part can not be observed because the 
flux is limited. Therefor, according to the PC method, 
the density at the center of binl is estimated as 10.55 
(the total source number in right part of binl divided 
by its four-volume) which is significantly smaller than 
the input model density. It can be expected that the 
more intensely the LF evolves with redshift, the more 
significant the error of (j>pc for binl is. The above dis- 
cussion also applies to the 1/Va method. 

2.7 A simple rule of thumb to determine the redshift 
and luminosity intervals 

The accuracy of binned LF also depends on how the 
redshift and luminosity bins are divided. The bins are 
supposed to be small enough, but can not be too small 
to include only a few objects. In general, the equal in- 
tervals of luminosity are used, e.g. AL ~ 0.3 or 0.5. 
While in the literature, the redshift intervals are com- 
monly chosen something arbitrarily, e.g. 0.0 < z < 0.2, 
0.2 < z < 0.5, 0.5 < z < 1.0; 0.0 < z < 0.15, 
0.15 < z < 0.4, 0.4 < z < 0.7 etc. We believe this 
may lead the bins located at the faint end to enclose 
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Fig. 10 Determination of the redshift and luminosity bins, 
(a): The redshift intervals are chosen arbitrarily, (b): The 
redshift intervals are determined by the intersecting points 
of luminosity grid-line and flux limit curve, (c): The lu- 
minosity intervals are determined by the intersecting points 
of redshift grid-line and flux limit curve. The shaded re- 
gions represent the surveyed regions. The surveyed regions 
only occupy very small parts of the faint end bins in (a), 
this probably lead them to enclose very few objects which 
gather at the right margins. This issue is avoided in (b) and 
(c). 

very few objects inside(see Fig. 10 (a)) and cause bias 
with small number statistics. Moreover, in these bins, 
objects necessarily gather at the right margins (see Fig. 
4(c)). Thus the area equal to [ALVq] of these bins will 
be much larger than that of the surveyed regions, caus- 
ing a significantly small estimate of 4>i/Va- To tackle 
this issue, we take the redshift intervals to be deter- 
mined by the intersecting points of luminosity grid-line 
and flux limit curve (see Fig. 10 (b)). Alternatively, if 
one persists in dividing redshift intervals randomly, the 
luminosity intervals should be determined by the inter- 
secting points of redshift grid-line and flux limit curve 
(see Fig. 10 (c)). 

Fig. 11 shows the binned LFs produced by the 1/Va 
and PC methods for our simulated data with the red- 
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Fig. 11 Binned luminosity functions of simulated samples of objects using (pi/Va 't'PC with the binning chosen according 
to the rule of thumb suggested in section 2.7. The input model LFs take values of z = 0.12,0.35,0.75,1.25,1.75,2.25 
respectively and are shown as dashed lines. 



shift and luminosity intervals determined according to 
the rule suggested above (Fig. 10 (c) specifically). It 
is noticed that adopting the rule of thumb indeed im- 
proved the binned LFs to some extent, especially for 

is mo re sensitive to h o w the 
Page fc Carreral (|200d ) 



opc- 



4>i/v^- Obviously, 
intervals are chosen than 
noticed that for objects close to the flux limit ipi/Va 
nearly always to be too small. We believe this is due to 
the arbitrary choosing of redshift and luminosity inter- 
vals. If the intervals are chosen appropriately, (t>i/Va 
only slightly larger than <j>pc at the faint end of LFs. 
And it is hard to say the estimate of PC method is 
markedly better than the 1/Va method. In this sense, 
the improvement of the PC method over the 1/Va 
method is probably slight. 



3 Non-parametric Methods 

A non-parametric method is necessary when investigat- 
ing the LFs, as it constructs the LFs directly from the 
data and makes no assumptions about the form of the 
LFs. Examples of non-para metric methods in clude the 
classical ^ = N/V me thod |ChristenseJl975[) . the 0/$ 
method (lTurner]|l979l). the Step - Wise Max Likelihood 



method ( Efstathiou et al.l I1988D and the C method 



( Lvnden-Bel]ll97ll ). as well as the 1/Va and PC meth- 
ods discussed here. 



The C~ method (also see IWillmeij 119971 : iFan et al 
2001 ; ChoIoniewskilll98"7 ) is believed to have the advan- 



tage over the 1/Va method as it does not require any 
binning of dat a, and rests on a. strong mathematical 
foundation fsee lWoodroofe [[19851 ). The key assumption 
in the C~ method is that the luminosities and redshifts 
are independent, but researches on the LFs of AGNs 
have indicated that there arc both luminosity and den- 
sity evolutions and the i ndependence assumption is 



incorrec t fe.g.. IPe] Eggst IWillott et al.l l200l[). Some 



autho rs (jEfron fc Petrosianll992HMalonev fc Pet: 



rosian 



19991 ) developed the C method to remove the corre- 



lation of luminosities and redshifts by defining new in- 
dependent variables (say L' = L/gk{z) and z, where 
the function gkiz) describes the luminosity evolution). 
Then the C^ method can be used for the (L', z) data. 
In this sense, the improved C~ method is no longer a 
non-parametric method as gkiz) is parameter depen- 
dent. Usually, a simple form gk{z) = (1 -I- z)*^ is used, 
where A: is a free parameter. Then a so-called test statis- 
tic r is performed to determine fc by making T(k) = 



(see 



Malonev fc Petrosianlll999HSingal et al.H201ll) . In 



most c ases, the (l-t-z)*^ form is too simplistic. E.g. JPei 
( 199,^ used a Gaussian form and two free parameters 
were introduced. Consequently, the unknown luminos- 
ity evolution gk{z) restricts the C~ method. 

Notably, in recent years more rigorous approaches, 
although not all of which are non-parametric, have been 
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proposed. E.g., ISchafeij (|2007( ) developed a powerful 
semi-parametric approach which is built on a nonpara- 
metric extension of maximum likelihood called local 
likelihood modeling. It is expected that a new method 
needs time to be recognized and in widespread use. Be- 
sides, when combining multiple samples to estimate the 
LFs, the truncation boundary (as discussed in iScliaferi 
20071) of data in the L~ z plane is further complicated. 
The new non-parametric methods are also supposed to 
deal with this challenge. 



4 Summary 

The classical 1/T4 and PC methods of constructing 
binned luminosity functions (LFs) are revisited and 
compared by graphical analysis. The two methods give 
different estimate of (f) for the bins which are crossed by 
the flux limit curves L = Lumiz). Us i ng the combined 
sample established by Yuan fc Wan j ( 2012 ). we show 
that (f>i/Va trends to give a smaller estimate than (j)pc 
at the faint end of LFs, while it gives a larger estimate 
in the middle of the LFs. 

Using a combined sample simulated by a Monte 
Carlo method, the estimate of two methods are com- 
pared with the input model LFs. The two methods 
give identical and ideal estimate for the high luminos- 
ity points of each rcdshift interval. However, for the low 
luminosity bins of all the redshift intervals both meth- 
ods give smaller estimate than the input model. We 
conclude that once the LF is evolving with redshift, 
the classical binned methods will unlikely give an ideal 
es timate over the total lu minosity range. 



Page fc Carreral (|2000[ ) noticed that for objects close 



to the flux limit (/fi/y^ nearly always to be too small. We 
believe this is due to the arbitrary choosing of redshift 
and luminosity intervals. We noticed that (pi/Va, is more 
sensitive to how the binning are chosen than (j)pc- We 
suggest a new binning method, which can improve the 
LFs produced by the 1/T4 method significantly, and 
also improve the LFs produced by the PC methods. 
Our simulations show that after adopting this new bin- 
ning, both the 1/Va and PC methods have comparable 
results. 
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